Magnetoelastic effects in doped Fe2P 
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We use combine high resolution neutron diffraction (HRPD) with density functional theory (DFT) 
to investigate the exchange striction mechanism at the Curie temperature {Tc) of Fe2P and to 
examine the effect of boron and carbon doping on the P site. We find a significant contraction of 
the basal plane on heating through Tc with a simultaneous increase of the c-axis that results in 
a small overall volume change of ^ 0.01%. At the magnetic transition the Fej-Fe/ distance drops 
significantly and becomes shorter than Fe/-Fe//. The shortest metal-metalloid (Fe/-P/) distance 
also decreases sharply. Our DFT model reveals the importance of the latter as this structural change 
causes a redistribution of the Fe/ moment along the c-axis (Fe-P chain). We are able to understand 
the site preference of the dopants, the effect of which can be linked to the increased moment on the 
Fe/-site, triggered by both valence electron and structural contributions. 

PACS numbers: 75.30.Sg,75.30.Kz, 75.80.+q, 75.30.Et 
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I. INTRODUCTION 



Fe2P-based magnetic alloys attract interest from many 
areas of physics, materials science and geophysics re- 
search. They are fomid in meteorites and are considered 
as a candidate minor phase present in the Earth's coreii^. 
An understanding of the mechanism of formation of these 
minerals can help to identify the histories of planetary 
bodies and the composition of the Earth's outer core. 
In LiFeP04-based battery materials, a percolating nano- 
network of metal-rich phosphides including Fe2P can sig- 
nificantly enhance electrical conductivity^. Fe2P-based 
alloys can furthermore be prepared as 1-dimensional (ID) 
nanowires and nanocables'*. Of most relevance to this ar- 
ticle, however, is the prospect of hexagonal Fe2P-based 
alloys being used as room temperature magnetic refrig- 
erants. 

Fe2P exhibits a first order magnetic transition from a 
ferromagnetic (FM) state to a paramagnetic one (PM) 
at 217 K'"' accompanied by a significant change in the 
c/a-ratio of the hexagonal structure. The magnetisa- 
tion deviates from the Curie- Weiss law for temperatures 
up to 700 K and net magnetism can be observed above 
the Curie temperature in fields of only a few Tesla. In 
nanocable form, the magnetic transition temperature is 
shifted 10-50 K higher compared to the parent composi- 
tion because of strong strain and/or carbon doping^. The 
increase in the Curie temperature, Tc with only a small, 
partial replacement of phosphorus with other p-block ele- 
ments (B, Si or As) is remarkable. 10% replacement of P 
by B leads to ^120% change, while the same amount of 
Si and As substitution also results in a ^70% and ~60% 
increase (Fig.[T]) respectively, with a simultaneous change 



in the nature of the transition from first order to second 
order^"^. Such large changes in Tc are not restricted to 
doping by p-block elements. Partial replacement of Fe 
by Mn results in a significant increase of the saturation 
magnetisation, while the first order nature of the meta- 
magnetic transition is preserved up to and beyond room 
temperature and is tuned by magnetic field at a rate of 

The metamagnetism of Fe2P shares many features 
with the itinerant electron metamagnetism (lEM) of 
La(Fe,Si)i3^*. That material also has a PM to FM transi- 
tion, the temperature of which can be tuned from around 
190 K to well above 300 K. The high magnetisation state 
can be induced by a magnetic field above Tc, and the 
metamagnetic transition is tuned by field at a rate of 
around 4 KT~^. As in Fe2P the first or second order 
nature of its metamagnetism depends on the substituent 
(e.g. Co or Mn) or intercolating atom (e.g. H or C) ^^'*"^ . 
Fe2P-based and La(Fe,Si)i3-based alloys are two of the 
leading contenders for scale- up as magnetic refrigerants, 
due to their tunable metamagnetism, their large room 
temperature magnetocaloric effect (MCE) and the fact 
that they are mostly composed of abundantly available 
3d and p-block elements^'*. 

There has been a recent growth in theoretical inves- 
tigations of the origin and tuneability of metamagne- 
tim in Fe2P and of the appearance of a body centered 
orthorhombic structure in substituted alloys. This is 
partly motivated by the sensitivity to doping of the Curie 
temperature and any associated thermomagnetic hystere- 
sis, and given added impetus more recently by the in- 
vestigation of the MCE of industrially-scaled quantities 
of material where good compositional tolerance is re- 
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quired. A Landau-Ginzburg free energy analysis based 
on fixed- spin- moment (FSM) calculations that took into 
account the effect of spin fluctuations revealed the meta- 
magnetic nature of the Fe-atoms at a particular crys- 
tallographic (3/) site in the parent alloy as well as in 
the doped counterpartsi^^. Using linear muffin tin or- 
bitals in the atomic sphere approximation (ASA), Severin 
et alji^ found that the magnetic moments in the inter- 
related hexagonal and orthorhombic structures are very 
similar and scale with the nearest-neighbour Fe-P dis- 
tances. We recently found that the phonon vibrational 
free energy stabilizes the hexagonal phase, whereas the 
electronic and magnetic entropies favour a low symme- 
try orthorhombic structure in Si-substituted Fe2Pi-KSix 
alloysAS. An analysis of the exchange constants in the 
hexagonal structure of B-, Si- and As-doped Fe2P helps 
to explain trends in Curie temperature. A principal in- 
terlayer Fe-Fe interaction was identified as controlling the 
strength of ferromagnetism^^. 

In the present paper we employ a joint experimental- 
theoretical approach that we have used elsewhere to un- 
cover the microscopic mechanism of metamagnetism in 
Mn-based antiferromagnets. By using density functional 
theory (DFT) in combination with structural informa- 
tion from high resolution neutron diffraction (HRPD)2£ 
we have previously mapped the magnetic phase diagram 
of Mn-based metallic orthorhombic compounds in the 
Pnma space group as a function of Mn-Mn distance^ and 
have predicted new antiferromagnetic metamagnets2Si2^. 
Taking the same approach here, we first analyse and com- 
pare the peculiar magnetoelastic coupling of the metam- 
agnetic transition in Fe2P, in boron-doped Fe2(P,B) and, 
for the first time, carbon-doped Fe2(P,C) using HRPD. 
In the second part of the article, we use a simple DFT 
model to further interpret the large interatomic changes 
at the magnetic transition. Our aim is to provide clues as 
to how to control and optimise the structure-magnetism 
relationship in this highly tuneable materials class. 



II. PREVIOUS EXPERIMENTAL FINDINGS 

Fe2P is the prototype structure of the hexagonal space 
group 189 (P62m) with 9 atoms in the unit cell. The 
6 Fe atoms occupy two non-equivalent threefold symme- 
try sites (3f and 3g), while the phosphor atoms sit on 
a singlefold (lb) and on a twofold position (2c) in the 
crystal lattice. Here we adopt the following established 
notation: Fe/ for the 3f (x/, 0, 0) positions, Fe// for the 
3g {xii, 0, 1/2) positions, and P/ and P// for the 2c 
(1/3, 2/3, 0) and lb (0,0, 1/2) positions of the P atoms, 
respectively. The hexagonal cell is composed of triangles 
in the ab plane as shown in Figure [21 The iron atoms 
that alternate along the c-axis are surrounded either by 
four P- atoms with tetrahedral symmetry (Fe/) or by 
five P-atoms forming a pyramid (Fe//). It was reported 
previously21 that the c-axis lattice parameter increases 
with temperature while the a- and 6-axes (basal plane) 




6.10- 
6,05 - 

< a 6.00- 

0) 5,95 - 
S 5 90 - 
(U 




X (at%) 



Figure 1: (Color online) The effect of boron, silicon and ar- 
senic doping in Fe2Pi-x'Zx on the magnetic ordering temper- 
ature Tc (top) and the room temperature lattice parameters 
(bottom)^"!. 




Figure 2: (Color online) Atomic arrangement of the Fe2P in 
the basal plane. Fe atoms occupy two non-equivalent three- 
fold symmetry sites, 3f (Fe/) and 3g (Fe//), and the phospho- 
rus atoms sit on a single- fold lb (P/) and two- fold position, 
2c (Pit) in the hexagonal crystal lattice. Fe/ {xi, 0, 0) atoms 
share the same plane with P/ (1/3, 2/3, 0) and Fe// atoms 
{xii, 0, 1/2) with P// (0,0, 1/2), respectively. 



exhibit negative thermal expansion. The same authors 
used a strain gauge dilatometer to measure the linear 
thermal expansion (Al/l) of a single crystal of Fe2P upon 
cooling through the first order type magnetic transition 
and found a sharp increase of = 0.74 x 10"'^ with the 
simultaneous decrease of the c-axis: (-^ — — 0.84xl0~^). 

The coupling between the structure and magnetism is 
further manifested by the effect of doping shown in Fig.[T] 
A significant increase in Curie temperature upon partial 
replacement of the P atoms by other p-block elements 
such as B, Si and As has been observed, Catalano et 
al.^ investigated the effect of isovalent substitution of As 
for P on the crystalline lattice and on magnetic proper- 
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ties. The authors found that the larger atomic radius 
As can be continuously accommodated into the hexago- 
nal lattice between < a; < 0.65 in Fe2Pi-a:As3; with a 
strong site preference for the P/ (2c) position. The mono- 
tonic increase of the lattice volume with As addition is a 
compromise of a moderate increase of a-axis and a simul- 
taneous contraction of the lattice along the c-direction. 
The same tendency was observed by the substitution of 
larger atomic radius, non-isovalent Si for P by Jernberg 
et al4 Interestingly, Si atoms were also found to show 
some site preference for the P/ position with opposing 
trends in the a and c lattice parameters and a signifi- 
cant increase in the magnetic ordering temperature. Sil- 
icon addition above x < 0.1 results in a change from the 
hexagonal lattice structure into one with a body-centered 
orthorhombic, Imm2 (44) symmetry. 

Finally, the partial replacement of P by the substan- 
tially smaller atomic radius boron results in a formidable 
increase of the Curie temperature. The solid solubility 
of the B atoms in Fe2Pi-a:B2; is, however, limited to 
X ~ 0.15 due to the formation of other Fe5PB2, Fe^B 
and Fe2B refractory borides at higher B concentrations. 
Chandra et al. established by means of Mossbauer spec- 
troscopy that the small boron atoms, unlike the larger 
elements (As and Si), occupy the P// (lb) singlefold po- 
sition in the hexagonal lattice^. Another consequence 
of the chemical pressure on the hexagonal lattice caused 
by boron addition as compared to the parent alloy or 
to the Si/ As-doped compositions is decreased lattice vol- 
ume. Despite such differences, the monotonic increase of 
the a and decrease of the c lattice parameters of the unit 
cell is strikingly similar to the effect of larger p-element 
dopants (Fig. [1]). An additional consequence of boron 
doping is that the first-order PM-FM transition of un- 
doped Fe2P becomes second order. 



III. METHODS 

A. Experimental Methods 

The samples used in this study were prepared from ul- 
tra high purity elements. The powders were mixed to 
weight ratios to according to nominal compositions of 
Fe2P and Fe2P0.96Z0.04 , (Z=B or C) and then thoroughly 
ground together in an agate mortar under protective at- 
mosphere. The initial powders were then pressed into 
pellets and sealed into quartz ampoules under protective 
argon atmosphere for solid state reaction. The initial 
annealing temperature was raised slowly (0.5 K min~^) 
up to 673 K, where each sample was kept for 4 hours 
for an initial reaction. There was a subsequent heating 
step (at 1 K min^^) up to 1273 K where the temperature 
was held for a further 4 hours. An additional heat treat- 
ment at 973 K for 4 hours was applied before the sample 
was oven-cooled to room temperature. X-ray diffractom- 
etry (XRD) was used to evaluate the structural prop- 
erties of the prepared samples. Single phase hexagonal 



compositions were found in all specimens by this method. 
Neutron diffraction was carried out at the time-of-flight 
high resolution powder diffractometer (HRPD) at ISIS, 
UK. This instrument has a resolution of ^ of 10^^ and 
was used at temperatures between 4.2 and 550 K. Neu- 
tron diffraction found some traces of unreacted carbon in 
Fe2P0.96C0.04 which we here refer to as a nominal compo- 
sition since the actual carbon concentration of the main 
phase will be somewhat smaller. Magnetic properties 
were measured between 10 and 400 K in a Quantum 
Design Physical Properties Measurement System (QD- 
PPMS). 



B. Computational model 



The electronic structure calculations were performed 
using the Vienna ab initio simulation package (VASP) 
code, based on DFT within projector augmented wave 
(PAW) method^s with Perdew-Burke-Ernzerhof (PBE) 
parameterization^^. Site-based magnetic moments were 
calculated using the Vosko-Wilk-Nusair interpolation^! 
within the general gradient approximation (GGA) for the 
exchange-correlation potential. A fc-point grid of 11 x 11 
X 13 was used to discretize the first Brillouin zone and 
the energy convergence criterion was set 10^^ eV during 
the energy minimization process. The density of states 
(DOS) plots presented in this work were calculated on a 
dense (19 x 19 x 21) fc-grid for high accuracy. The spin- 
orbit interaction was turned off during the calculations 
and only coUinear ferromagnetic (FM) and non-magnetic 
(NM) configurations were considered. 

The minimal, nine-atom basis cell (six Fe atoms and 
three P atoms) was used to evaluate the total energies 
and magnetic properties of the alloys. Using this sim- 
ple model, the effect of doping was simulated by the 
replacement of a single phosphorous atom by another 
p-block element (Z) that represents an x=l/3 compo- 
sitional change in the Fe2Pi_a;Zj; formula. Although this 
approach is undoubtedly oversimplified in respect of ex- 
act compositions provided in the experimental section, 
we believe it is still a suitable model to capture the rel- 
evant changes in the electronic structure caused by the 
doping elements. In order to be consistent with the ex- 
perimental results, we only considered changes along the 
a- and c-axis by the individual expansion and compres- 
sion of the a- and c-lattice parameters, without allow- 
ing any relaxation of the strained structure. In prac- 
tice, we varied the lattice parameters using a stepsize 
of ±0.5%, calculating the self-consistent electronic struc- 
tures at each step. 
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Figure 3: (Color online) Thermomagnetic curves of the parent 
Fe2P alloy with the C and B doped counterparts (top, left 
axis) in the field of ^oAf=0.01T. Derivatives (^^) are linked 
to the right axis (top). Bottom figure shows the change of 
the lattice parameters (a, c) for comparision, obtained using 
HRPD. 
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Figure 4: (Color online) Temperature evolution of the lattice 
parameters of the parent Fe2P compound together with that 
of the C- and B-doped samples. The strong magnetoelastic 
response is especially apparent around the magnetic ordering 
temperature (~215 K). 



B. High resolution neutron diffraction 



IV. RESULTS 

A. Magnetometry 

The samples were first cooled to 10 K, where mag- 
netisation curves were cohected at fields up to 5 T, after 
which thermo-magnetic curves were collected on heat- 
ing to 400 K in an applied field of 0.01 T. The results 
are plotted in Fig. \3\ The parent alloy shows a sharp, 
first order transition at T sa215 K in accordance with 
values reported previously. The strong effect of doping 
on the magnetic order temperature is apparent. As ex- 
pected, the partial replacement of P by the much smaller 
atomic radius B significantly increases Tq- In this study, 
we also partially replaced P by C atoms for the first 
time. The empirical atomic radius of carbon is much 
smaller (70 pm) than that of the phosphorus (100 pm) 
and boron (85 pm) elements. On the other hand, in terms 
of valence electrons the sequence B(p^) <C{p^) <P{p^) 
stands. In practice, carbon doping shows a very sim- 
ilar effect to that of the other p-block substituents as 
it also clearly increases the magnetic order tempera- 
ture. Furthermore, the Curie transition is heavily broad- 
ened by doping as demonstrated by the smearing out of 
the temperature derivative of the magnetisation {^^) 
(Fig. [3]). Finally, the saturation magnetisation {Ms) in- 
creases slightly with doping; Ms=112 , 113 and 116 ^ 
was obtained for Fe2P, Fe2P0.96C0.04 and Fe2P0.96B0.04j 
respectivelly at 10 K in 9 T applied magnetic field. This 
slight increase of Ms is in line with the expectations from 
DFT calculations (see sec. IV Bl) . However, the samples 
were not fully saturated due to the large magnetocrys- 
talline anisotropj^. 



Fig. m shows the anisotropic lattice expansion of the 
a- and c-axes of Fe2P and the doped compounds. In 
the magnetically ordered state (T <200 K) the basal 
plane exhibits negative lattice expansion with increas- 
ing temperature. When the temperature reaches Tc, a 
sharp contraction of the lattice in the ab plane is ob- 
served. The a-lattice expansion only looks Debye-like at 
higher temperatures (T >240 K) in the paramagnetic 
phase. On the other hand, the thermal expansion of the 
c-axis is found to be positive over the entire investigated 
temperature range. The magnetic ordering temperature 
is also strongly reflected in the lattice response along c. 
The consequence of these counteracting lattice parameter 
changes over the magnetic transition is a volume change 
at Tc which is as small as 0.01%. 

The effect of doping on the lattice parameters is clear; 
tlie basal plane of the hexagonal lattice expands whilst 
tlie c-axis shows contraction. At first glance, the anoma- 
lous thermal lattice expansion resembles that of the par- 
ent alloy. However, temperature derivatives of these 
quantities {da/dT, dc/dT) reveal characteristic differ- 
ences between the parent and doped alloys in Fig. [3l The 
sharpness of the derivatives demonstrates the first order 
nature of the magnetoelastic transition of stoichiometric 
Fe2P. Although da/dT and dc/dT are of opposite sign, 
they both have sharp peaks at temperature Tp=209 K. 
Both B and C doping cause a shift of Tp to higher tem- 
peratures with broader da/dT and dc/dT. These effects 
on the lattice properties correspond very well to the ob- 
served changes in the magnetic properties. The temper- 
ature derivative of the magnetisation in Fig. [3] (top), re- 
veals similar broadening of the ferromagnetic transition 
with B or C doping. 

HRPD is capable of tracking the change in the inter- 
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Figure 5: (Color online) The shortest Fei-2Fe/ and Fei- 
2Fe// distances as a function of temperature for all inves- 
tigated samples. Inset depicts the temperature region of the 
crossover. 



atomic distances, thus providing vital information about 
the magnetoelastic coupling in these materials. The evo- 
lution of the lattice parameters through the magnetoe- 
lastic transition is already a clear indication of the sig- 
nificant changes of the atomic distances. Both the 3f and 
3g positions of the iron atoms are low symmetry posi- 
tions described by the positional parameters xi and xu, 
respectively. Iron atoms on the 3f-site (Fe/) are con- 
nected to 2 other iron atoms of the same type, denoted 
as Fe/-2Fe/. There are also two distinctive Fe/-2Fe// 
and Fe/-4Fe7j connections to iron atoms located at the 
3g-site (Fe//). Finally, there exists a relatively close Fe//- 
4Fe// distance above 3 A. In Fig. [SI we only plot the two 
shortest distances, Fe/-2Fe/ and Fe/-2Fe//, as a function 
of temperature, evaluated from Rietveld refinements of 
the HRPD data. In all of the investigated samples, the 
Fe/-2Fe// distance is the shortest Fe-Fe separation at low 
temperatures(<200 K). The latter increases with tem- 
perature, eventually becoming larger than the Fe/-Fe/ 
distance, which is strongly reduced in the vicinity of the 
magnetic transition temperature. The point in tempera- 
ture above which Fe/-Fe/ is the shortest Fe-Fe distance 
is between 190 and 220 K for all samples. 

In addition we can distinguish two groups of metal- 
metalloid distances: Fe/-(P/,P//) and Fe//-(P/,P//). 
The shortest one is Fe/-2P/, followed by Fe/-2P//, Fe//- 
P// and Fe//-4P/. Fig. [H] contains the two shortest dis- 
tances only. The separation of Fe//-P// and Fe//-4P/ 
atoms is in the range of ^2.37 A and ^2.49 A, respec- 
tively. It is worth noting that both Fe/-2P/ and Fe//- 
P// distances only have components within the a6-plane 
which explains the resemblance of their temperature evo- 
lution to that of the a lattice parameter in Fig. 0) the 
shortest Fe/-2P/ distance decreases sharply around the 
magnetic transition temperature in the parent alloy and 



< 

<u 
o 
sz 
m 

(/> 

Q. 
0) 



2.305 

2.300 

2.295. 

2.290 

2.285 

2.280: 
2.220. 

2.218 

2.216' 

2.214. 

2.212 

2.210 




■Fe P B 

2 0.96 0.04 

■Fe P C 

2 0.96 0.04 



Fe-P, 



50 100 150 200 250 300 350 400 450 

Temperature (K) 

Figure 6: (Color online) The shortest metal-metalloid dis- 
tances, Fe/- 2P/ and Fe/- 2P// as a function of temperature 
in Fe2P together with the doped counterparts. Both distances 
relate to the metamagnetic 3f (Fe/) site. 
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Figure 7: Total magnetisation as a function of lattice param- 
eter of Fe2P. Magnetisation is strongly linked to the change 
in a-lattice parameter. 



the doped compounds. The larger a lattice parameter 
in the doped compounds means that the metal-metalloid 
distance is largest in those samples. 



V. THEORETICAL RESULTS 

A. Stoicheometric Fe2P 

The variation of total magnetisation with lattice pa- 
rameters expected from our DFT calculations is plotted 
in Figure[7l A change in the lattice can have a very strong 
effect on the magnetic properties. The lattice expansion 
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Figure 8: (Color online) Comparision of the densitiy of states 
of a high magnetisation state (red) and a low magnetisation 
state (black) in Fe2P. 



within the basal plane causes the total magnetisation to 
increase abruptly from ^ 1.3 /is to 3.1 fiB around a crit- 
ical value of a w 5.7 A, and stays practically constant 
above it. The magnetisation varies very little with c, re- 
flecting the greater importance of the atomic distances 
within the basal plane, where the atoms are packed more 
densely. In order to depict the changes related to the 
electronic structure that are most relevant to such an 
anisotropic magneto-elastic coupling, we compare the to- 
tal electronic density of states (DOS) of the high and low 
magnetisation states, at (a, c) values indicated by red and 
black dots in Figure [71 The lattice structure in the high 
magnetic state has a high and sharp DOS at Ep with 
'-^80% spin polarisation, dominated by minority spins 
(Figure [51 red line). This signature is typical of a mag- 
netic instablility. When the lattice is compressed within 
the basal plane, the magnetic ground state is significantly 
different, as shown in black in Figure [SI The filled bands 
in the bottom contain the phosphorus 3s states. The 
conduction band is formed by the 3d states of Fe and the 
3p states of P as is typical of strong p-d hybridization. 
The hybridized states are also present near the Fermi 
level, influencing the exchange splitting of the 3d states 
of Fe. The effect of the smaller a-lattice parameter (low 
magnetisation state) is to reduce the exchange splitting 
considerably. The high peaks in the majority density of 
states located around -0.5 eV are shifted to -0.2 eV. In 
the minority DOS, the opposite trend is observed as the 
triple peak feature between and 0.25 eV is lowered to 
around -0.25 eV. 

The relationship between the DOS at the Fermi En- 
ergy, D{Ep) and the strength of magnetic exchange, /, 
can be expressed as the Stoner criteria / x D{Ef) > 1 
for an itinerant ferromagnet. Such an analysis has been 
extensively discussed by Eriksson et alj^*. They inves- 



tigated site-related magnetic stability and also found a 
magnetic instability with decreasing cell volume. In our 
study we concentrate on comparing the differences in spin 
densities of the high and low magnetic configurations, 
plotted in Fig. [21 As a consequence of the lattice symme- 
try, the shortest metal- metalloid distance, Fej-P//, as 
well as the shortest metal- metal separation, Fe/-Fe//, 
have projections along both the a- and c-axes, and so 
this is the most suitable cross-section for our analysis. 

We can ascribe three competing effects to the mag- 
netoelastic transition. Firstly, with increasing tempera- 
ture the separation between Fe/ and Fe// atoms increases 
monotonously as observed experimentally (see Fig. [SJ. 
This is linked to our recent investigations of the exchange 
coupling mechanism among the magnetic sites using the 
exact muffin-tin orbital (EMTO) method. Our exchange 
coupling analysisiS. found that the largest contribution 
to the decrease in the c/a ratio to the total magnetic 
exchange interaction is the weakened magnetic interac- 
tion between the Fe atoms on 3f and 3g sites. Secondly, 
at the magnetic transition the Fe/-Fe/ distance is also al- 
tered and it becomes comparably to the Fe/ and Fe// but 
the exchange constant between these atoms is only the 
fraction of the Fe/ and Fe// exchange, and has a less pro- 
nounced dependence on distance. Thirdly, the shortest 
Fe/-P/ distance, located in the ab plane also dips around 
the transition temperature. This results in a significantly 
stronger hybridisation that can explain the decrease in 
magnetisation and alternation of these bonds, and hence 
trigger a lattice distortion in the first place. A similar 
scenario was also drawn from the comparison of the total 
electronic charge density in the FM and PM states of Mn 
and Si doped Fe2P-type alloy by Dung et ali^ 

Our current DFT results provide further evidence for 
an understanding of the microscopies of the Curie tran- 
sition. The calculations reveal the unusual duality of 
the magnetic structure of Fe2P in which there is a 
large magnetic moment of M(Fe//)=2.16 /is on the 3g- 
site together that coexists with a significantly smaller, 
Af(Fe/)=0.85 /is, moment on the 3f crystallographic site. 
These findings are in full agreement with previous stud- 
ies (see Sec. [lT|. If we consider solely the effect of the 
structural changes that we know to occur at the Curie 
transition (Fig. [5]) , we find that there is a redistribution 
of magnetisation of the Fe/-site. The large and localised 
magnetic moment on the Fe// site is decreased in ampli- 
tude to 1.3/LtB above the transition and strong delocal- 
isation of the Fe/ magnetisation occurs simultaneously. 
Electrons from the Fe/-site that have a magnetic con- 
tribution are redistributed along the c-axis in directions 
that link the Fe/ and P// sites. This strong coupling be- 
tween magnetism and lattice is the origin of the strong 
magnetoelastic response and the highly field dependent 
transition in these alloys. In the following section, we 
compare the effect of doping on the structure and mag- 
netic properties obtained by means of DFT using struc- 
tural data from neutron diffraction (HRPD). 
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Figure 9: (Color online) A comparison of the spin density of 
the high magnetisation state (right) and low magnetisation 
state (left) of Fe2P. The large magnetic moment around the 
Fe/j site is only decreased in amplitude above the transition, 
when strong delocalisation of the Fej magnetisation occurs 
simultaneously. 



B. The Effect of Doping 

There exist two crystallographic sites where the metal- 
loid elements can be incorporated into the lattice, albeit 
with strong site preference as described in Sec.|lll Exper- 
imentally, elements with an atomic radius smaller than P 
were found to show a preference in the single-fold position 
(lb), whilst larger elements tends to occupy the twofold 
symmetry site (2c) in the hexagonal lattice. In order to 
establish the site preference from DFT calculations, we 
compare the difference in total energy for the different 
elemental substitution. First of all, the calculated total 
energies of the FM solutions (at any (a, c) lattice param- 
eter value) were always found to be energetically more 
stable than that of the NM configurations by typically 
1.5-3 eV/f.u.; yielding a strong tendency toward the for- 
mation of magnetic order in all the investigated alloys. 
The site preference of the dopants can be established by 
calculating the difference in total energy {dE) between 
the single-fold and two-fold site position occupancies us- 
ing the equilibrium lattice of the FM state. The calcula- 
tions find a clear trend with the size of the substituents: 
the energy difference is negative for elements that are 
smaller than P (left hand side of Fig. [TOl) , reflecting the 
preference for singlefold occupation. On the other hand, 
a substituent with larger atomic radii (Si and As) prefers 
to occupy the twofold position of the hexagonal lattice. 

This theoretically established site occupation of the 
dopant elements is in line with the experimental find- 
ings see Secini suggesting the dominance of size effects. 
Therefore, we hereafter focus our magnetoelastic investi- 
gations on the effect of single-fold site occupation by the 
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Figure 10: (Color online) Site preference of the p-block 
dopants as established from DFT. Negative total energy dif- 
ferences represent preference for the single-fold position. Pos- 
itive total energy differences suggest tendency for two-fold site 
occupancy. 
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Figure 11: (Color online) Equilibrium lattice parameters for 
NM and FM states as a result of DFT calculations as a func- 
tion of dopants in Fe2P2/3Zi/3, where Z=C, B, Si, or As, 
respectively. 



small C and B elements and the two-fold site occupation 
of the large Si and As elements. 

The effect of doping on the lattice parameters is shown 
in Fig. [TlJ The non-magnetic (NM) calculations (black 
line) show an increasing trend in both the a- and c pa- 
rameters of the hexagonal structure with the inclusion of 
larger atomic radius, anionic p-block elements. This be- 
haviour is easily anticipated in terms of chemical pressure 
as the lattice is expected to adapt according to the size 
of dopant. Examining Fig. [TTl the steepest change in 
the lattice is seen with the partial replacement of P by B 
atoms, due to the largest difference in atom size between 
these elements. However, in the presence of ferromag- 
netic (FM) interactions, any monotonic relation of lattice 
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Z= C B P Si As 



Figure 12: (Color online) Total and site-projected magnetic 
moments in ¥62? 2/3'^!/ 31 for Z=C,B,Si and As, respectively. 



parameter to the atomic radius of the dopant is broken. 
The parent alloy has significantly smaller a and larger c 
parameters in the FM ground state than the doped com- 
pounds. The onset of ferromagnetism increases the sepa- 
ration of atoms within the basal plane (positive exchange 
striction) while their separation along the c-axis is re- 
duced (negative exchange striction) , regardless of the size 
or valence state of the dopant as compared to the parent 
composition (Fe2P). The overall volume change caused 
by the ferromagnetic exchange is positive, dominated by 
the larger a parameter {Vhex = a^csin(27r/3)). Based 
on these calculations, one would expect the lattice to ex- 
pand along the a-axis with a simultaneous contraction of 
the c-axis in the vicinity of the magnetic ordering tem- 
perature. Indeed, experimental (HRPD) measurements 
clearly reveal this sharp positive exchange striction in the 
basal plane and the shrinkage of the c-axis upon cooling 
the sample through the Curie temperature (see Fig.|3]). 

The significant coupling of the lattice dimensions to 
the magnetism also becomes aparent from the analy- 
sis of the total and site-projected magnetic moments 
as plotted in Fig. [T^l The magnetic moment of Fe// 
atoms show very little variation with doping and stays 
at around 2.2 hb- Although the larger p-block dopants 
increase the separation of the Fe//-Fe// atoms, this dis- 
tance is typically above 3 A, and so doping only slightly 
alters the exchange-split states. In strong contrast, the 
site-resolved magnetisation of the 3f site is found to be 
much more sensitive to the valence electron number of 
the anionic elements. Examining Fig. [1^ there is a 
strong relation between added valence electron number 
and magnetisation with Z=B (— 2e~), C (— le^), P(— ), 
Si (— 2e^), or As(— ). This feature can be explained in 
terms of charge transfer from Fe/ to the anionic Z element 
in order to fill the electronegative p shell of the latter. 
The charge transfer is distance dependent and is strongly 
supressed by the strong p-d interaction at smaller Fe/-Z 
separations, while lattice expansion suppresses p-d hy- 



bridization and increases the degree of localisation and 
ionic bonding in the system. With doping, the short 
Fe/-Z distances are strongly influenced, altering the bal- 
ance between magnetic moment formation (Fig. [5]) and 
bond formation. The larger magnetic moments on the 
Fe/ atoms that result from doping account for the expan- 
sion of the basal plane as indicated by the comparison of 
Figs, [m and El 



VI. SUMMARY AND CONCLUSION 

We have used high resolution neutron diffraction 
(HRPD) as well as density functional theory to investi- 
gate the effect of p-block elements doping on the magne- 
toelastic properties of Fe2P. HRPD has revealed that the 
strong coupling between the magnetism and the lattice is 
manifested by the contraction of the basal plane and by 
a significant increase of the c-axis on heating through the 
magnetic transition, resulting in a small overall volume 
change of the lattice (>0.01%). A simultaneous change 
in both the metal-metal and metal-metalloid distances is 
observed. 

DFT calculations reveal that the magnetic properties 
are strongly dependent on expansion in the basal (ab) 
plane, while they are almost invariant with regard to 
variation of the c parameter. The strong dependence 
of magnetic moments on the 3f site is related to the a 
lattice parameter being on the verge of a metamagnetic 
transition. The closest metal-metal and metal-metalloid 
distances - the latter ones also linked to the metamag- 
netic 3f site - are strongly altered by both the d-d and 
p-d hybridisation energies at the transition. As a result, 
the redistribution (delocalisation) of the magnetisation 
from the 3f site along the Fe/-P// chains in the c-axis 
direction occurs as shown in Fig. ^ implying its strong 
influence on bonding (also suggested by Dung et a\.— ). 
The magnetic moment on the 3g site is only slightly re- 
duced (to ~ 1.3/zb) above the transition, accounting for 
the observed strong divergence of the Curie- Weiss law^. 

Our theoretical investigations on the effect of doping 
show that elements with an atomic radius smaller than 
phosphorus (e.g. C or B) occupy single-fold (lb) sites 
while larger elements such as Si or As prefer to occupy 
the two-fold symmetry site (2c) in the hexagonal lattice. 
The magnetism on the high magnetic site (3g) is not in- 
fluenced by the dopant and stays approximately constant 
(^2.2/iB). By contrast, the iron atoms on the metam- 
agnetic site develop larger magnetic moments, in direct 
relation to the valence electron number of the doping el- 
ement. As a result, the conditions for the metamagnetic 
transition of the 3f site are locally modified, resulting 
in a smearing of the transition as observed experimen- 
tally. Indeed, the larger exchange splitting on Fe/ results 
in increased exchange coupling parameters, significantly 
increasing the Curie temperature^-. 

The large difference in the lattice equilibrium between 
the non-magnetic and ferromagnetic solutions in the par- 
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ent compound as well as in the doped materials as found 
by DFT calculations suggests that the magnetic proper- 
ties are not only very sensitive to compositional inhomo- 
geneities but also to internal strains set by the prepa- 
ration conditions'^. Therefore for applications (such as 
magnetic refrigeration), where the transition tempera- 
ture needs to be set by the material to a high accuracy, 
elements such as Si and As are favorable but process con- 
trol will be crucial. Finally, as the remarkable change in 
magnetic ordering of these alloys is linked to the lattice 
parameters within the basal plane, uniaxial thin films on 
flexible substrates could be exploited to electrically ma- 
nipulate the magnetic properties of this material system. 
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